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Abstract 

We construct the Liouville operator for the SU(2) classical colored Coulomb plasma (cQGP) for 
arbitrary values of the Coulomb coupling T = V/K, the ratio of the mean Coulomb to kinetic en- 
ergy. We show that its resolvent in the classical colored phase space obeys a hierarchy of equations. 
We use a free streaming approximation to close the hierarchy and derive an integral equation for 
the time-dependent structure factor. Its reduction by projection yields hydrodynamical equations 
in the long- wavelength limit. We discuss the character of the hydrodynamical modes at strong 
coupling. The shear viscosity is shown to exhibit a minimum at T ~ 8 near the liquid point. This 
minimum follows from the cross-over between the single particle collisional regime which drops as 
1/r 5 / 2 and the hydrodynamical collisional regime which rises as T 1 / 2 . The self-diffusion constant 
drops as l/r 3//2 irrespective of the regime. We compare our results to molecular dynamics simu- 
lations of the SU(2) colored Coulomb plasma. We also discuss the relevance of our results for the 
quantum and strongly coupled quark gluon plasma (sQGP). 
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I. INTRODUCTION 



High temperature QCD is expected to asymptote a weakly coupled Coulomb plasma 
albeit with still strong infrared divergences. The latters cause its magnetic sector to be non- 
perturbative at all temperatures. At intermediate temperatures of relevance to heavy-ion 
collider experiments, the electric sector is believed to be strongly coupled. 

Recently, Shuryak and Zahed |]J have suggested that certain aspects of the quak-gluon 
plasma in range of temperatures (1 — 3) T c can be understood by a stronger Coulomb in- 
teraction causing persistent correlations in singlet and colored channels. As a result the 
quark and gluon plasma is more a liquid than a gas at intermediate temperatures. A liquid 
plasma should exhibit shorter mean-free paths and stronger color dissipation, both of which 
are supported by the current experiments at RHIC ^]. 

To help understand transport and dissipation in the strongly coupled quark gluon plasma, 
a classical model of the colored plasma was suggested in [3j. The model consists of mas- 
sive quarks and gluons interacting via classical colored Coulomb interactions. The color 

lowing from Poisson brackets. For the 



is assumed classical with all equations of motion fo 

SU(2) version both molecular dynamics simulations fl and bu.a thermodynamics flQ were 
recently presented including simulations of the energy loss of heavy quarks jg]. 

In this paper we extend our recent equilibrium analysis of the static properties of the 
colored Coulomb plasma, to transport. In section 2 we discuss the classical equations of 
motion in the SU(2) colored phase space and derive the pertinent Liouville operator. In 
section 3, we show that the resolvent of the Liouville operator obeys a hierarchy of equations 
in the SU(2) phase space. In section 4 we derive an integral equation for the time-dependent 
structure factor by introducing a non-local self-energy kernel in phase space. In section 5, we 
close the Liouville hierarchy through a free streaming approximation on the 4-point resolvent 
and derive the self-energy kernel in closed form. In section 6, we project the self-energy 
kernel and the non-static structure factor onto the colorless hydrodynamical phase space. 
In section 7, we show that the sound and plasmon mode are the leading hydrodynamical 
modes in the SU(2) colored Coulomb plasma. In section we analyze the shear viscosity for 
the transverse sound mode for arbitrary values of T. We show that a minimum forms at 
r w 5 at the cross-over between the hydrodynamical and single-particle regimes. In section 
8, we analyze self-diffusion in phase space, and derive an explicit expression for the diffusion 
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constant at strong coupling. Our conclusions and prospects are in section 9. In appendix A 
we briefly summarize our variables in the SU(2) phase space. In appendix B we detail the 
projection method for the self-energy kernel used in the text. In appendix C we show that 
the collisional color contribution to the Liouville operator drops in the self-energy kernel. In 
appendix D some useful aspects of the hydrodynamical projection method are outlined. 



II. COLORED LIOUVILLE OPERATOR 

The canonical approach to the colored Coulomb plasma was discussed in [3| . In brief, the 
Hamiltonian for a single species of constituent quarks or gluons in the SU(2) representation 
is defined as 

N 2 N n n 

i 2m * 4ii r <- p ii 

The charge g 2 /4:7i has been omitted for simplicity of the notation flow and will be reinserted 
in the pertinent physical quantities by inspection. 

The equations of motion in phase space follows from the classical Poisson brackets. In 
particular 



^ = -{H,r i } = ^Lp. = ?i (11.2) 
dt dpj ovj m 

The Newtonian equation of motion is just the colored electric Lorentz force 



*~W*>~gg-flT*-F. ("3) 
with the colored electric field and potentials defined as (a = 1,2,3) 

O a 

El = -V^ = -V< > , 3 , (II.4) 

Our strongly coupled colored plasma is mostly electric following the original assumptions 
m . The equation of motion of the color charges is 



M = _ WQ . } = _ E g^ { 05 iQ ; } = E ^ („,, 

j,b J j^i ■" 
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for arbitrary color representation. For SU(2) the classical color charge (111.51) precesses around 
the net colored potential $ determined by the other particles as defined in (1II4I) . 



^ = (** x Qd (II.6) 



This equation was initially derived by Wong [8|. Some aspects of the SU(2) phase space are 
briefly recalled in Appendix A. 

The set ( 111.21) . ( 1II.3|) and ( 111.51) define the canonical evolution in phase space. The time- 
dependent phase distribution is formally given by 



A? 



/(*, rpQ) = £ *(r - nOWip ~ Pi(t))S(Q - QM) = E *(« - ( IL? ) 



i=l 



For simplicity q is generic for r,p, Q. Using the chain rule, the time-evolution operator on 
(pX7]) obeys 

d _ d dvj d dpj d dQ i d _ 

dt dt + dt dr t + dt dp, + dt dQ, -° t + lL - ^> 
The last relation defines the Liouville operator 

C = C + d + C Q = -iv i -V ri -i Fi ■ V Pi - ■ (Qi x V Qi ) (II.9) 

The last contribution in (III. 9[) is genuily a 3-body force because of the cross product (orbital 
color operator). It requires 3 distinct colors to not vanish. This observation will be important 
in simplifying the color dynamics below. Also (111.91) is hermitean. 
Since (111.7j) depends implicitly on time, we can write formally 



j t f(t,r P Q)=ijCf(t,rpQ) (11.10) 

with a solution f(t) = e lCt f(0). The formal relation (IH.lOj) should be considered with 
care since the action of the Liouville operator on the 1-body phase space distribution (1II.7j) 



generates also a 2-body phase space distribution. Indeed, while Co is local in phase space 
£o £ 5(q - g 4 ) = -iv ■ V r £ % - g 4 ) = Lofo) £ ~ ft) ( IL11 ) 
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the 2 other contributions are not. Specifically 



i=Fj,mn J 

- J dq'Lj(q, q') ]T % - 9m ) W - qrj (11.12) 



with 



L z (9, q') = iVr|^^| ■ (V p - VpO (11.13) 



Similarly 



E ^ _ «m) = ~ l E ^ ~ ^, | ■ V Q« ^ " 9m) 

= - /" ^(^E^-gJ^g'-gJ ( IL14 ) 



with 



£g(g, g') = j^T77| ■ ( v q - V Q') ( IL15 ) 

Clearly fill. 141) drops from 2-body and symmetric phase space distributions. It does not for 
3-body and higher. 



III. LIOUVILLE HIERARCHY 

An important correlation function in the analysis of the colored Coulomb plasma is the 
time dependent structure factor or 2-body correlation in the color phase space 



S (t - t ',r- r', pp', Q Q') = (5f(t, rpQ) 5f(t', r'p'Q')) 
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(III.l) 



with Sf = f — (f) the shifted 1-body phase space distribution. The averaging in (IHI.lj) 



is carried over the initial conditions with fixed number of particles iV and average energy 
or temperature (3 = 1/T. Thus (/) = nfo(p) which is the Maxwellian distribution for 
constituent quarks or gluons. In equilibrium, the averaging in (fillip is time and space 
translational invariant as well as color rotational invariant. 
Using the ket notation with 1 = q = rpQ 



IW,i)> 

with also 2 = q', 3 = q", 4 = q'" and so on and the formal Liouville solution 5f(t, 1) = 
e lCt 8f{l) we can write (IHI.ip as 

S(t-t',q,q') = (5f(t,l)|5f(t',2)> = (5f(l)\e iC ^\6f(2)) (III.3) 

The bra-ket notation is short for the initial or equilibrium average. Its Laplace or causal 
transform reads 

S(z,q,q') = -i / dt9(t-t')e lzt S(t-t',q,q') = (5f{l)\——\6f(2)) (III.4) 

J-oo Z + L 

with z = uj + iO. Clearly 



zS(z,q,q') + ( < y/(l)| J C_L-| ( y/(2)> = (5/(l)|5/(2)) (III.5) 

Z + £ 



Since £t = C is hermitian and using (jll.lip . (III. 1 2[) and (III. 14[) it follows that 



(5/(l)|£ = (5f(l)\L (q) - J dq» L I+Q (q,q") (5/(l)5/(3)| (III.6) 



Thus 



(z-Loiqtys&qrf)- J dq" 'L I+Q (q> \q")S(z,qq" \q') = S (q 7 q') (LIU) 
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where we have defined the 3-body phase space resolvent 



S(z,qq",q') = (5/(1)5/(3) | 



z + C 



\Sf(2) > 



(III.8) 



Sq(QiQ') is t ne static colored structure factor discussed by us in [9j. Since Lj + q(q' , q) is 
odd under the switch q <-> q', and since S(z, qq", q') = S(—z, q, q'q") owing to the t <-> t' 
in pILijl . then 



(z + L Q {q'))${z,q,q')- J dq"L I+Q (q', q")S(z, q, q'q") = S (q,q') (111.9) 



( \lll.7\i or equivalently (1III.9P define the Liouville hierarchy, whereby the 2-body phase space 
distribution ties to the 3-body phase space distribution and so on. Indeed, 1 )11191) for instance 
implies 



z + 



L (q"))$(z, qq', q") - J dq"'L I+Q {q", q"')S(z, qq', q"q"') = S (qq', q") (111.10) 



with the 4-point resolvent function 



S(z,qq',q"q"')= (5f(l)5f(2)\- 



1 



15/(3)5/(4) > 



(iii.li; 



[ z + C 

These are the microscopic kinetic equations for the color phase space distributions. They 
are only useful when closed, that is by a truncation as we discuss below. These formal 



equations where initially discussed in 
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Il2l . Il3j in the context of the one component 



Coulomb Abelian Coulomb plasma. We have now generalized them to the multi-component 
and non- Abelian colored Coulomb plasma. 

IV. SELF-ENERGY KERNEL 



In (IIII.7j) the non-local part of the Liouville operator plays the role of a non-local self- 
energy kernel S on the 2-body resolvent. Indeed, we can rewrite (1III.7j) as 



[z - L (q))S(z, q, q') - J dq"Z(z, q, q")S(z, q" , q') = S (q, q') 



(IV. 1) 
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with the non-local self-energy kernel defined formally as 



J dq"Y>(z, q, q") S(z, q" , q') = J dq"L I+Q (q, q") S(z, qq", q') (IV.2) 

The self-energy kernel £ can be regarded as the sum of a static or z-independent contribution 
Ss ans a non-static or collisional contribution Sc, 

E(z, q, q") = E s (q, q") + X c (z, q, q") (IV.3) 
The stationary part £5 satisfies 

J dq"E s {q, q") S (q", q') = J dq" L I+Q {q, q") S (q, q' , q") {IV A) 

which identifies it with the sum of the 2- and 3-body part of the Liouville operator Lj + q. 

The collisional part is more involved. To unwind it, we operate with (z + L (q')) on 
both sides of flIV.2j) . and then reduce the left hand side contribution using (IIII.9I) and the 
right hand side contribution using (1III.10I) . The outcome reduces to 

E c (z, q, q") 8 (q", q') = -J dq'" L I+Q (q, q")L I+Q {q\ q") S(z, qq" , q'q'") 

+ J dq'" E(z, q, q") L I+Q (q', q") S(z, q" , q'q"') (IV.5) 

after using flIV.41) . From flIV.21) it follows formally that 

Q") = J dq'" L I+Q (q, q'") S" 1 ^, q\ q") S(z, qq'" , q') (IV.6) 

Inserting (1IV.6[) into the right hand side of (1IV.5[) and taking the q' integration on both sides 
yield 

nf (p") E c (z, q, q") = - J dq" dq"' L I+Q (q, q")L I+Q (q', q'") G(z, qq" , q'q'") (IV. 7) 



S 



with G a 4-point phase space correlation function 



Giz^q'^q'q'^ =S(z,qq[,q'q' 2 )- J dq 3 dq 4 S(z, qq[, q 3 )S \z, q 3 , g 4 )S(z, q 4 , q'q 2 ) (IV.8) 



The collisional character of the self-energy Sc is manifest in (jIV.71) . The formal relation for 



the collisional self-energy (1IV.7I) was initially derived in 12|, ll3[] for the one-component and 
Abelian Coulomb plasma. We now have shown that it holds for any non-Abelian SU(N) 
Coulomb plasma. 

Eq. (1IV.7I) shows that the connected part of the self-energy kernel is actually tied to a 4- 
point correlator in the colored phase space. In terms of flIV.7j) . the original kinetic equation 
(IIII.7j) now reads 



(z - L (q)) 8{z, q, q') - j dq"j: s (q, q")S(z, q\ q') = S (q, q') 

- J dq" dq x dq 2 L I+Q (q, q^Lj+Qiq", q 2 ) G(z, qq^ q"q 2 ) S(z, q" , q') (IV.9) 

which is a Boltzman-like equation. The key difference is that it involves correlation functions 
and the Boltzman-like kernel in the right-hand side is not a scattering amplitude but rather 
a reduced 4-point correlation function. ( 1IV.9I) reduces to the Boltzman equation for weak 
coupling. An alternative derivation of ( 1IV.9I) can be found in Appendix C through a direct 
projection of flIV.2p in phase space. 



V. FREE STREAMING APPROXIMATION 

The formal kinetic equation ( 1IV.7I) can be closed by approximating the 4-point correlation 
function in the color phase space by a product of 2-point correlation function [3] , 



G(t, <?<h, q'q 2 ) « (S(t, q, q')S{t, q u q 2 ) + S(t, q, q 2 )S{t, q' , (V.l) 

This reduction will be referred to as the free steaming approximation. Next we substitue 
the colored Coulomb potentials in the double Liouville operator Li+q x L\ + q with a bare 
Coulomb V(r -r',Q Q') = Q Q'/\r - r'\. 
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Li+q(Q, Qi) = *V r V(r -r x ,Q Q x ) • (V p - V Pl ) 

-i (Q x V Q V(r -r^Q-QJ-VQ + Q.x V Ql V(r - n, Q QJ V Ql ) (V.2) 
times a dressed colored Coulomb potential defined in J9] 



1. 



L? +Q (q, q x ) = -i^V r c D (r - n, Q • (V p - V Pl , 



(Q x V Q c D (r - n, Q Q x ) ■ V Q + Q x x V Ql c D (r - n, Q • V Ql ) (V.3) 



14j in the context of the 



This bare-dressed or half renormalization was initially suggested 
one-component Coulomb plasma to overcome the shortcomings of a full or dressed-dressed 
renormalization initially suggested in 12|, ll3). The latter was shown to upset the initial 
conditions. Thus 



Li + Q(q,qi)L I+Q (q',q 2 ) -»■ - (L I+Q (q, q 1 )Lf +Q (q', q 2 ) + Lf +Q (q, q x )L I+Q {q' , q 2 )^ (V.4) 



Combining (jVT]l and (jV^jl in pVTFj) yields 



" /o(p') Sc(*, 9, 9') « -- y dqfi dq 2 ( At+q(<Z, qi)Lf +Q (q', q 2 ) S (^ 9; <z') S (*> 9i, 9 a ) 
qi)Lf +Q (q', q 2 )S(t, q, q 2 )S(t, q\ q x ) + (<?i ^ q 2 ,q^ q')) (V.5) 

This is the half dressed but free streaming approximation for the connected part of the 
self-energy for the colored Coulomb plasma. Translational invariance in space and rota- 
tional invariance in color space allows a further reduction of (IV. 5 j) by Fourier and Legendre 
transforms respectively. Indeed, Eq. (1V.5|) yields 
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~ y dflfi dflfa \Li{q,qi)Lf(q',q 2 )S(t,q, q')S{t, q x ,q 2 ) 
+Li(q, q x )Lf{q', q 2 )S(t, q, q 2 )S(t, q', q x ) + «-> q 2 , q q')^j 

= ~^J d <li d <l2 (v r c D (r - n, Q Q x ) V p V r ,V(r' - r 2 , Q' ■ Q 2 ) ■ V^S(t, q, g')S(*, 9i, 9 2 ) 

+V r c D (r - n, Q ■ V p V r ,V(r' - r 2 , Q' ■ Q 2 ) • V p >S(t, g, g 2 )S(t, <?', + ^q 2 ,q^ q'U 

(V.6) 

where we note that the colored part of the Liouville operator dropped from the collision 
kernel in the free streaming approximation as we detail in Appendix C. Both sides of flB.6j) 
can be now Legendre transformed in color to give 



n/o(f/) ^Zci&rr'^rt^—PtiQ.Q') 

i 

» -± ^ti(i±Lp I+1 (Q . Q<) + ^fl-W • Q')) 

x ( V r c D1 (r - ri) • V p V r /— , • V^S^t, rr', pp') s i(^ r ir 2 , P1P2) 

\ \r' - r 2 \ 

+V r c Dl (r - n) ■ V p V r /— f • V p /Si(t,rr2,pp2) s «(^ r ' r i ; P'Pi) 

r - r 2 \ 

+V r ,c m (r' - r 2 ) ■ V p 'V rl — - — , • V p Si{t,rr 2 ,pp 2 )S 1 (t,r'r 1 ,p'p 1 ) 

\r — r*i| 

V r /c D1 (r' -r 2 ) • V p ,V r ^^ • V p S,(t, rr', pp'^t, r ir2 ,p lP2 ) J (V.7) 



Thus 
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nf (p') E C i{t,rr',pp') 

~ -77^ / dr 1 dp 1 dr 2 dp 2 ( V r c Di (r - n) ■ V p V r /— r • Vjy 



2/ + 

+V r c D i(r - ri) ■ V p V r /— r ■ V p / 



x (^/"Vy S '- l( ^' PP')Si(*, rir 2 , PxPa) + 2/+T S,+1 ^' rr ' pp ') Sl ^' rir2 ' P1P2)) 



F - 1*2 1 

x(^PYS 1 (t,rr2,pp 2 )Si-i(t,r , V 1 ,p'p 1 ) + ^^-Si(t,rr a ,pp 2 )Si + i(t,rVi,p'p 1 ) 

/ / 1 
+V r /c D ;(r - r 2 ) • V p 'V rl r ■ V p 

\r — ri I 

x (^PYS z _i(t, rr 2) pp 2 ) Si (^rV^p'pi) + ^^-S/ + i(t,rr 2 ,pp 2 )Si(t,rVi ) p'pi) 
+V r 'C m (r' - r 2 ) ■ VwV rl — - — r • V p 

1 7* — 7*i| 

X (27TI S ^ l( ^' rr " pp '^ Sl ^ rir2 ' PlP ^ + ^T Sz+1 ^' rr '' pp '^ Sl ^^ riT " 2 ' PiPa))l 

(V.8) 

with S^_i = by definition. In the colored Coulomb plasma the collisional contributions 
diagonalize in the color projected channels labelled by I, with I = being the density channel, 
/ = 1 the plasmon channel and so on. In momentum space (1V.8I) reads 



nfo(p') Zci(t,k,pp' 

1 r . . r di 



2(3 



J dp±dp 2 J (l ■ V p l ■ Vj/c D1 {l)Vi 



x (2TTi S ^ 1 ^' k ~ l ' pp "> Sl ^ Z 'PiP2) + 2^T S ' +1 ^' k ~ ^pp')Si(t, ^,pip 2 ) 

+l-V p (k-l)-V p ,c Dl (l)V k _ l 

x (2T+^ Sl<yt, k ~ Z 'ff2)Si-i(t,i,p / Pi) + ^-j 81 ^' k ~ ^pp 2 )Si+i(t,z,p / Pi) 

+(k - I) • Vpl • Vp,c Dl (l)V k -i 

X (^~j^ S '- 1 ^' Z ' pp 2)SK t > k - hp'Pi) + ^yS i+1 (t,Z,pp 2 )Si(t, fc - l,p'p 1 )^j 
+l-V p l-Vp>c D1 {l)Vi 

x (2TTT S ^ 1 ^' fc ~~ ^pp')Si(t, z,pip 2 ) + ^r^+ifo k - *,pp')Si(t, z,pip 2 ))^ 

(V.9) 
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with Vi = An /I 2 . We note that for I = which is the colorless density channel (IV.9P involves 
only Si which is the time-dependent charged form factor due to the Coulomb interactions. 



VI. HYDRODYNAMICAL PROJECTION 



In terms of fTV9l) , flTV2l and 



EiOzfc.ppj.) = ( Eqi + Sjj + S Q i + Ecj (zk,pp 1 



(VI. 1) 



the Fourier and Legendre transform of the kinetic equation flIII.7j) now read 



zS t (zk,pp') - J dpJl^zk^pp^S^zk^p-tp') = S ol (k,pp') 



(VI.2) 



with S ; = L and Y> S i = L(i+q)i- Specifically 



Ejj^fe.ppi) = -n/ (p) c Di (fc) 

m 

E ol (^fe,p Pl ) = (VI.3) 

and Ec; is defined in ( IV.9[) . See also Appendix B for an alternative but equivalent derivation 
using the operator projection method. 

(1VI.2j) is the key kinetic equation for the colored Coulomb plasma. It still contains 
considerable information in phase space. A special limit of the classical phase space is the 
long wavelength or hydrodynamical limit. In this limit, only few moments of the phase 
space fluctuations 5f or equivalently their correlations in S ~ (Sfdf) will be of interest. In 
particular, 



n(t, r) 
p(t,r) 
e(t,r) 



dpdQ 5f(t,r,p, Q) 
dpdQp5f(t,r,p,Q) 



dpdQ^-6f(t,r,p,Q) 
2m 



(VI.4) 
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The local particle density, 3- momentum and energy (kinetic). The hydrodynamical sector 
described by the macro-variables flVI.4j) is colorless. An interesting macro-variable which 
carries charge representation of SU(2) would be 



m(t,r) = / drdQYr(Q)5f(t,r,p,Q) (VI.5) 

m 

which reduces to the I color density with I = being the particle density, I = 1 the charged 
color monopole density, I = 2 the charged color quadrupole density and so on. Because of 
color rotational invariance in the SU(2) colored Coulomb plasma, the constitutive equations 
for (1VI.5j) which amount to charge conservation hold for each /. 

To project flVI.2j) onto the hydrodynamical part of the phase space characterized by flVI.5j) 



and (1VI.4I) . we define the hydrodynamical projectors 



Pi? = X}|i><i| Qh = 1h-V h (VI.6) 

i=i 

with 1 = 1-density, 2,4,5 = momentum and 3 = energy as detailed in Appendix D. When 
the / = particle density is retained in ( 1VI.6I) the projection is on the colorless sector of the 
phase space. When the I = 1 charged monopole density is retained in flVI.6j) the projection is 



on the plasmon channel, and so on. Most of the discussion to follow will focus on projecting 
on the canonical hydrodynamical phase space (IVI.4I) with I = or singlet representation. 
The inclusion of the I ^ representations of SU(2) is straightforward. 
Formally ( 1VI.1I) can be viewed as a p x p 1 matrix in momentum space 



(z-Etizk)) S l (zk) = S m (k) (VI.7) 

The projection of the matrix equation (1VI.7j) follows the same procedure as in Appendix B. 
The result is 

(z - Vu^i{zk)V H - V H Qi{zk)V H ) V H Si(zk)V H = V H S ol (k)V H (VI.8) 

with 

6, = mzk)Q H (z - QH^Hiz^QH^Qn^zk) (VI.9) 
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If we define the hydrodynamical matrix elements 



G Uj (zk) = (i\Si(zk)\j) 
Z Uj (zk) = (i|£,(*fc)|j) 
e Uj (zk) = (i\&i(zk)\j) 

G olij (zk) = (i\S ol (k)\j) (VI.10) 



then NTM reads 



(z5 u > - Quj(zk)) Giji^zk) = Goiu'(fc) 



(VI.11) 



with Qi = E; + ©«• (IVI-lip takes the form of a dispersion for each color partial wave / with 
the projection operator (lVI.6j) set by the pertinent density (lVI.5j) . The contribution to 
Qi will be referred to as direct while the contribution will be referred to as indirect. 



VII. HYDRODYNAMICAL MODES 

The zeros of ( 1VI.11I) are the hydrodynamical modes originating from the Liouville equa- 
tion for the time-dependent structure factor. The equation is closed under the free streaming 
approximation with half renormalized vertices as we detailed above. 

We start by analyzing the 2 transverse modes with i = T in fIVI.lOj) and (I VI. lip . We note 
with [15|] that Gin = whenever T ^ i. The hydrodynamical projection (see Appendix D) 
causes the integrand to be odd whatever I. The 2 independent transverse modes in fIVI.lip 
decouple from the longitudinal i = L, the (kinetic) energy i = E and particle density i = N 
modes for all color projections. Thus 



Gi T (zk) = — — - (VII.l) 

z — \Li T {zk) 

with fli T = (T\Qi\T) and G/ T = (T\Gi\T). The hydro-projected time-dependent / structure 
factor for fixed frequency z = u> + iO, wavenumber k develops 2 transverse poles 



zi{k) = Qi T (zk) w 0{k 2 ) 
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(VII.2) 



The last estimate follows from 0(3) momentum symmetry under statistical averaging what- 
ever the color projection. We identify the transverse poles in (I VI 1. 21) with 2 shear modes of 
consititutive dispersion 



uo + i-^k 2 + 0(k 3 ) = (VII.3) 
mn 

with t]i the shear viscosity for the 1th color representation. Unlike conventional plasmas, the 
classical SU(2) color Coulomb plasma admits an infinite hierarchy of shear modes for each 
representation I. 




(c) (d) 



FIG. 1: Soi(q) from SU(2) Molecular Dynamics. 
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The remaining 3 hydrodynamical modes L,E,N are more involved as they mix in (IVI.llj) 
and under general symmetry consideration. Indeed current conservation, ties the L mode 
to the N mode for instance. Most of the symmetry arguments regarding the generic nature 
of Qi in [15( carry to our case for each color representation. Thus, for the 3 remaining 
non-transverse modes (IVLllft reads in matrix form 



/ 



V 



GiNN 


GiNL 


G 


GiLN 


Gill 


G 


GlEN 


Giel 


G 



INE 
ILE 
IEE 



\ 



\ 



z 











—&ILE 





— &IEL 


Z ~ &IEE 



\ 



J 



7i 



v 



-nhi 
1 
1 



\ 



The 3 remaining hydrodynamical modes are the zeros of the determinant 



(VII.4) 



A, 



z -Vti NL {zk) 

-Qi LN (zk) z - Q lLL (zk) -Vt lLE {zk) 
-QiEL(zk) z - Q lEE (zk) 



(VII.5) 



( 1VII.5I) admits infinitly many solutions zi(k). We seek the hydrodynamical solutions as 
analytical solutions in k for small k, ie. z(k) = zi n k n for each SU(2) color representation 
I. In leading order, we have 



A z « zio z l0 



k 2 T 



m 



(VII.6) 



after using the symmetry properties of Qi as in [15( for each I. We have also made use of 



the generalized Ornstein-Zernicke equations for each I representation [9| 
the molecular dynamics simulation results for 4 typical structure factors 



n Fig. [T] we show 

3 



47T 



21 + 1 



E 



jrn 



3 ifc-£C,(0) ym 



Y l m (Q l 



(VII.7) 



for I = 0, 1, 2, 3. We have made use of the dimensionless wavenumber q = kaws with aws 
is the Wigner-size radius. In Fig. [2] we show the analytical result for Sqi which we will use 
for the numerical estimates below. We note that the I = 1 structure factor which amounts 
to the monopole structure factor vanishes at k = 0. All other Vs are finite at k = with 
I = corresponding to the density structure factor. 
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(1VII.6j) displays 3 hydrodynamical zeros as k — > for each / representation. One is mass- 
less and we identify it with the diffusive heat mode. The molecular dynamics simulations 
of the structure factors in Fig. [T] implies that all I ^ channels are sound dominated with 
two massless modes, while the / = 1 is plasmon dominated with two massive longitudinal 
plasmon states. Thus 

z l± = ±u 2 p 8 n (VII.8) 

with uj p = ko ^T/m the plasmon frequency. The relevance of this channel to the energy 
loss has been discussed in [16|. We used S i(fc ~ 0) ~ k 2 /k 2 D with k 2 D the squared Debye 
momentum. All even I ^ 1 are contaminated by the sound modes. The SU(2) classical 
and colored Coulomb plasma supports plasmon oscillations even at strong coupling. These 
modes are important in the attenuation of soft monopole color oscilations. 



1.4- 




VIII. SHEAR VISCOSITY 

The transport parameters associated to the SU(2) classical and colored Coulomb plasma 
follows from the hydrodynamical projection and expansion discussed above. This includes, 

18 



the heat diffusion coefficient, the transverse shear viscosity and the longitudinal plasmon 
frequency and damping parameters. In this section, we discuss explicitly the shear viscosity 
coefficient for the SU(2) colored Coulomb plasma. 

Throughout, we define A = |7r(3r) 3//2 , the bare Coulomb interaction = k\jl 2 in units 
of the Wigner-size radius k$ l = aws- While varying the Coulomb coupling 



r = (VIII.l) 

4tt a W s 

all length scales will be measured in aws — (47rn/3) _1/ ' 3 , all times in the inverse plasmon 
frequency l/u p with uj 2 = K 2 D /m(3 = ng 2 C2/m. All units of mass will be measured in m. 
The Debye momentum is k 2 d = g 2 n(5Ci and the plasma density is n. for instance, the shear 
viscosity will be expressed in fixed dimensionless units of 77* = nmio v a ws . 

The transverse shear viscosity follows from (I VII. II) with contributing to the direct or 
hydrodynamical part, and B; contributing to the indirect or single-particle part. For I = 



V0_ _ Vodir VOind (VI 1 1 2) 

respectively. The direct or hydrodynamical contribution is likely to be dominant at strong 
coupling, while the indirect or single-particle contribution is likely to take over at weak 
coupling. We now proceed to show that. 

The indirect contribution to the viscosity follows from the contribution outside the hy- 
drodynamical subspace through Qh and lumps the single-particle phase contributions. It 
involves the inversion of Qh^coQh in (1B.13j) with 



mn |(t|S |t/)| 2 mn | (t|(S 00 + S co )|t/) | 2 

r/oind = hm - = lrm — — — — — (VIII. 3) 

fc^o k 2 (tl\iE \tl) k^o k 2 (tl\iT,co\tl) 

In short we expand T,co in terms of generalized Hermite polynomials, with the first term 
identified with the stress tensor due to the projection operator (ID. 31) . The inversion follows 
by means of the first Sonine polynomial expansion. Explicitly 



r] . d= Vo M = nm ]im l l(t|S 00 + S C0 (fc,0)|Zt)| 2 = (1 + A/ 2 ) 2 
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with 



1 1 



^ O °rfg(2(S 01 (g) 2 -l) + (l-S 01 (g))) 



60tt 2 (3ry/2 

J 3 = 777^77^ rdqq(l - S i(g)) (VIII.5) 



10vr 3 / 2 3r J 

with the dimensionless wave number q = kaws- 

We recall that Soi is the monopole structure factor discussed in j^] both analytically 
and numerically. In Fig. [2] we show the behavior of the static monopole structure factor 
from [9J for different Coulomb couplings. The larger F the stronger the first peak, and the 
oscillations. These features characterize the onset of the crystalline structure in the SU(2) 
colored Coulomb plasma. A good fit to Fig. [2] follows from the following parametrization 



1 + C e-^ Cl sin {{q - C 2 )/C 3 ) (VIII.6) 



with 4 parameters 6*0,1,2,3- The fit following from (1VIII.6I) extends to q ~ 100 within 10~ 5 
accuracy, thanks to the exponent. 

The direct contribution to the shear viscosity follows from similar arguments. From 
(I VI I. II) and (1VII.3I) . we have in the zero momentum limit 



TfLTl TflTl 

Vodir = lim — <t|i£ |t) = lim —(t\iZ co {0, 0)\t) (VIII.7) 

with S = Soo + S/o + Sco as defined in ( 1VI.3I) and ( IV. 91) . Only those nonvanishing contribu- 
tions after the hydrodynamical projection were retained in the second equalities in flVIII.3j) 
as we detail in Appendix D. A rerun of the arguments yields 



V* dir = Vo dir /v* = Jim ^jj^l dtn(e • if 

x (c m (i)g nl (k - z, t)g nl (i, t)v t - c m (i)g nl (k - z, t)g nl (i, t)v k ^ 



(VIII. 

The projected non-static structure factor is 
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g nl (i,t) = - f dpd P 's 1 (i,t;pp') = g nl (i,t)s 01 (i) (viii.9) 

n J 

with the normalization Q nl (l,0) = 1. As in the one component Coulomb plasma studied 
in [nj we will approximate the dynamical part by its intermediate time-behavior where the 
motion is free. This consists in solving (lIV.ip with no self-energy kernel or X = 0, 



g nl (l,t)^e-W 2 / 2m ?S 01 (l) (VIII.10) 

Thus inserting (IVIII.lOj ) and performing the integrations with k — > yield the direct con- 
tribution to the shear viscosity 



= — = zrla r * (vni.il) 

The full shear viscosity result is then 

Tjo = VQdir + Voind = p l + (1 + A/ 2 ) 2 j 

rj* rf rj* 457T 1 / 2 A/3 



after inserting flVIII.4jl and fIVIII.llh in (1VIII.2I) . The result (1VIII.12D for the shear viscosity 



of the transverse sound mode is analogous to the result for the sound velocity in the one 
component plasma derived initially in [14] with two differences: 1/ The SU(2) Casimir in T; 
2/ the occurrence of Soi instead of Soo- Since Soi is plasmon dominated at low momentum, 
we conclude that the shear viscosity is dominated by rescattering against the SU(2) plasmon 
modes in the cQGP. 

Using the fitted monopole structure factor (1VIII.6I) in ( 1VIII.5I) we can numerically assess 



( 1VIII.4P for different values of T. Combining this result for the indirect viscosity together 



with (IVIII.llI) for the direct viscosity yield the colorless or sound viscosity r/ . The values 
of 770 are displayed in Table I, and shown in Fig. [3] (black). The SU(2) molecular dynamics 
simulations in [3| which are parameterized as 



D 242 f) D72 

O.OOir + ^ + (VIII.13) 
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are also displayed in Table I and shown in Fig. [3] (red) for comparison. The sound viscosity 
dips at about r ~ 8 in our analytical estimate. To understand the origin of the minimum, 
we display in Fig. H] the scaling with T of the direct or hydrodynamical and the indirect 
part of the shear viscosity. The direct contribution to the viscosity grows like T 1//2 , the 
indirect contribution drops like 1/T 5 / 2 . The latter dominates at weak coupling, while the 
former dominates at strong coupling. This is indeed expected, since the direct part is the 
contribution from the hydrodynamical part of the phase space, while the indirect part is 
the contribution from the non-hydrodynamical or single-particle part of phase space. The 
crossing is at T rs 4. 
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FIG. 3: The direct and indirect part of the viscosity 



TABLE I: Reduced shear viscosity. See text. 

r 2 4 6 8 10 12 14 16 18 
rj* QGP 0.286 0.092 0.067 0.066 0.070 0.076 0.081 0.087 0.092 
T)md 0.217 0.168 0.168 0.139 0.132 0.127 0.124 0.122 0.120 



The reduced sound velocity 77* is dimensionless. To restore dimensionality and compare 
with expectations for an SU(2) colored Coulomb plasma, we first note that the particle 
density is about 3 x 0.244 T 3 = 0.732 T 3 . There are 3 physical gluons, each carrying 
black-body density. The corresponding Wigner-Seitz radius is then aws — (3/47™) 1 / 3 w 
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FIG. 4: The best fit of the direct and indirect part of the viscosity 

0.688/T. The Coulomb coupling is T ~ 1.453 {g 2 N c / At:). Since the plasmon frequency is 
u 2 = K 2 D /m(3 = ng 2 N c /m, we get u 2 ~ 3.066 T 2 (g 2 N c /Air) with m ~ 3T. The unit of 
viscosity t]q = nmcUpO^g translates to 1.822 T 3 (g 2 N c / Air) 1 / 2 . In these units, the viscosity 
for the SU(2) cQGP dips at about 0.066 which is r]* QGP « 0.066 r? « 0.120 T 3 (g 2 N c /A7r) 1 / 2 . 
Since the entropy in our case is o = 6 (47r 2 /90)T 3 , we have for the SU(2) ratio rj/a\su(2) = 
0.046 (g 2 N c / An) 1 / 2 . The minimum in the viscosity occurs at Y = 1.453 (g 2 N c /An) « 8, so 
that (g 2 N c /Air) 1 / 2 ~ 2.347. Thus, our shear viscosity to entropy ratio is rj/a\su{2) — 0.107. 
A rerun of these estimates for SU(3) yields r)/a\su(3) — 0.078 which is lower than the bound 
77/cr = 1/4-7T ~ 0.0795 suggested from holography. 

Finally, we show in Fig.[5]the shear viscosity T r j*s( q ) at low T (a:green) and large T (b:black) 
assessed using the weak-coupling structure factor S(k) = k 2 /(k 2 + k 2 D ). The discrepancy is 
noticeable for V near the liquid point. The large discrepancy for small values of V reflects 
on the fact that the integrals in (1VIII.5I) are infrared sensitive. The sensitivity is tamed by 
our analytical structure factor and the simulations. We recall that in weak coupling, the 
Landau viscosity t/l is [3] 



VL = I (VIII.14) 

77* 18 r 5 /2 \ n (r D /r ) 

which follows from a mean-field analysis of the kinetic equation with the plasma dielectric 
constant set to 1. The logarithmic dependence in (1VIII.14I ) reflects on the infrared and 
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FIG. 5: Comparison with weak coupling. See text. 

ultraviolet sensitivity of the mean-field approximation. Typically ro = 1/&d and vq 
(g 2 C 2 /4:7i)l3 which are the Debye length and the the distance of closest approach. Thus 



Vl 



5V37T 1 



(VIII.15) 



27 r 5 /2hi(i/r) 

or tjl/t]* rs 0.6/(r 5 / 2 ln(l/r)) which is overall consistent with our analysis. 

The Landau or mean-field result is smaller for the viscosity than the result from pertur- 
bative QCD. Indeed, the unsealed Landau viscosity (1VIII.15I) reads 



10 



Vl 



(VIII.16) 



nmiOpO^s and using ln(r D /r ) w 31n(l/a s )/2 with 



24 (a s C 2 ) 2 (3 5 / 2 a s 

after restoring the viscosity unit rj* 
a s = g 2 /4:TT. While our consituent gluons carry m ~ ttT, in the mean field or weak coupling 
we can set their masses to m ~ gT. With this in mind, and setting C 2 = N c = 3 in ( 1VIII.16I) 
we obtain 



Vl 



5V2 



T 3 



1087rV4 a 7/4 ln(1/ ^ 



0.05 



T 3 



al //4 ln(l/a s 



(VIII.17) 
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which is to be compared with the QCD weak coupling result [1 

Vqcd « 21 %r-, i (VIII.18) 
afln(l/a s ) 

The mean-field result ( 1VIII.17I) is ctj^ ~ yfg smaller in weak coupling than the QCD per- 
turbative result. The reason is the fact that in perturbative QCD the viscosity is not only 
caused by collisions with the underlying parton constituents, but also quantum recombina- 
tions and decays. These latter effects are absent in our classical QGP. 

IX. DIFFUSION CONSTANT 

The calculation of the diffusion constant in the SU(2) plasma is similar to that of the 
shear viscosity. The governing equation is again (1III.7P with X and S replaced by S s , S s . 
The label is short for single particle. The difference between S and S s is the substitution of 
(JED by 

f s (rpQt) = VN6(r - n(t))5(j> - Pl (t))6(Q - Q x (t)) (IX. 1) 

The diffusion constant follows from the velocity auto-correlator 



V D (t)= 1 -(V(t)-V(0)) (IX.2) 



through 



D= / dtV D (t) (IX.3) 



Solving (1III.7I) using the method of one-Sonine polynomial approximation as in Jj]] yields 
the Langevin-like equation 



(lh>lf) ^ dt'M(t')V D (t-t') (IX.4) 



o 



dt 

with the memory kernel tied to S^ , 
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nf (p') X s cl (t,k,pp') 

1 f f dl 

dp x dp 2 I j—^l • V p / ■ V p 'C m (/)Vj 



x 



(2Tn S ^ l(t ' fc ~ W)Si(M,PiP 2 ) + ^i; s f+i^ fe - hpp')Si(t,i, PlP2 ) 



(IX.5) 



and 



n/ (p') £g (t,fc = 0,pp') 
1 f f dl 

= --j d Pl dp 2 j —l ■ V P Z ■ Vp,c D1 (/)^Sf(t,/,pp , )Si(t,i,PiP 2 ) (IX.6) 



therefore 



M(0 = A /" rfpdp'p . p 'sg (t, fc = 0,pp')/o(p') (IX.7) 

which clearly projects out the singlet color contribution. If we introduce the dimensionless 
diffusion constant, D* = Djw v a^ vs , then (1IX.3I) together with ( 1IX.4I) yield 

If 00 1 r°° M(t) r°° 

— = m(3 dtMCt) -> — = 3r / w p dt^r^ = 3T / drM(r) (IX.8) 
D Jo D* J w 2 p Jo 

Using similar steps as for the derivation of the viscosity, we can unwind the self-energy kernel 

S s in flDCHD to give 



1 f dl f 00 

w = ~ v j w?L dri2 ° m ^ v ^ 1 ^ ^ 1 ^ (ix - 9) 

where we have used the same the half-renormalization method discussed above for the vis- 
cosity. The color integrations are done by Legendre transforms. Here again, we separate the 
time-dependent structure factors as Q n i(l,t) = S i(/){? n i(M) and Sj^^t) = Q n i(l,t) in the 
free particle approximation. Thus 



IF = r3 ' 2 (s) ' /°° dqq(l ~ Sol(,)) 
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(IX.10) 
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r 

FIG. 6: Diffusion Constant (black, green) versus molecular dynamics simulations (red). See text. 

TABLE II: Diffusion constant. See text. 

r 2 4 6 8 10 12 14 16 18 
D* QGP 0.410 0.115 0.055 0.034 0.024 0.017 0.014 0.012 0.010 
Di rn 0.230 0.132 0.095 0.076 0.063 0.055 0.048 0.044 0.040 



The results following from ( 1IX.10I) are displayed in Table [TT] and in Fig. [6] (black) from 



weak to strong coupling. For comparison, we also show the the diffusion constant measured 
using molecular dynamics simulations with an SU(2) colored Coulomb plasma The 
molecular dynamics simulations are fitted to 

For comparison, we also show the diffusion constant (lIX.lOp assessed using the weak coupling 
or Debye structure factor S(k) = k 2 /(k 2 + k 2 D ) in Fig. [H] (green). The discrepancy between 
the analytical results at small V are similar to the ones we noted above for the shear viscosity. 
In our correctly resummed structure factor of Fig. [2j the infrared behavior of the cQGP is 
controlled in contrast to the simple Debye structure factor. 
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Finally, a comparison of (lIX.lQj) to (1VIII.5j) shows that 1/D* ~ I/AJ3 which is seen to 
grow like T 3//2 . Thus D* drops like I/T 3 / 2 which is close to the numerically generated result 
fitted in Fig. [7J (left). The weak coupling self- diffusion coefficient scales as 1/r 5 / 2 as shown 
in Fig. [7J (right). More importantly, the diffusion constant in the SU(2) colored Coulomb 
plasma is caused solely by the non hydrodynamical modes or single particle collisions in our 
analysis. It does not survive at strong coupling where most of the losses are caused by the 
collective sound and/or plasmon modes. This result is in contrast with the shear viscosity 
we discussed above, where the hydrodynamical modes level it off at large T. 
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FIG. 7: Fit to the diffusion constant. See text. 



X. CONCLUSIONS 

We have provided a general framework for discussing non-perturbative many-body dy- 
namics in the colored SU(2) Coulomb plasma introduced in [l]. The framework extends the 
analysis developed intially for one-component Abelian plasmas to the non-Abelian case. In 
the latter, the Liouville operator is supplemented by a color precessing contribution that 
contributes to the connected part of the self-energy kernel. 

The many-body content of the SU(2) colored Coulomb plasma are best captured by 
the Liouville equation in phase space in the form of an eigenvalue-like equation. Standard 
projected perturbation theory like analysis around the static phase space distributions yield 
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a resummed self energy kernel in closed form. Translational space invariance and rigid color 
rotational invariance in phase space simplifies the nature of the kernel. 

In the hydrodynamical limit, the phase space projected equations for the time- dependent 
and resummed structure factor displays both transverse and longitudinal hydrodynamical 
modes. The shear viscosity and longitudinal diffusion constant are expressed explicitly in 
terms of the resummed self-energy kernel. The latter is directly tied with the interacting 
part of the Liouville operator in color space. We have shown that in the free streaming ap- 
proximation and half-renormalized Liouville operators, the transport parameters are finite. 

We have explicitly derived the shear viscosity and longitudinal diffusion constant of the 
SU(2) colored Coulomb plasma in terms of the monopole static structure factor and the for 
all values of the classical Coulomb parameter T — V/K, the ratio of the potential to kinetic 
energy per particle. The results compare fairly with molecular dynamics simulations for 
SU(2). 

The longitudinal diffusion constant is found to drop from weak to strong coupling like 
l/r 3//2 . The shear viscosity is found to reach a minimum for T of about 8. The large 
increase at weak coupling is the result of the large mean free paths and encoded in the 
direct or driving part of the connected self-energy. The minimum at intermediate T is tied 
with the onset of hydrodynamics which reflects on the liquid nature of the colored Coulomb 
plasma in this regime. 

At larger values of T an SU(2) crystal forms as reported in Our current analysis 
should be able to account for the emergence of elasticities, with in particular an elastic 
shear mode. This point will be pursued in a future investigation. The many body analysis 
presented in this work treats the color degrees of freedom as massive constituents with a 
finite mass and a classical SU(2) color charge. The dynamical analysis is fully non-classical. 
In a way, quantum mechanics is assumed to generate the constituent degrees of freedom 
with their assigned parameters. While this picture is supported by perturbation theory at 
very weak coupling, its justification at strong coupling is by no means established. 
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APPENDIX A: SU(2) COLOR PHASE SPACE 

A useful parametrization of the SU(2) color phase space is through the canonical variables 
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2l| 



Q 1 = COS 01 J 2 -71-2, Q2 = gin ( p i ^Jj2 _ 7r 2 ; q3 = ^1 ( Al ) 

with Q 2 being a constraint variable fixed by J 2 or the quadratic Casimir with q 2 = 
J2a C ~ 1 Q a Q a - The conjugate set Q l ,n l obeys standard Poisson bracket. The associated 
phase space measure is 



dQ = c R dnid(j)iJdJ5{J 2 - q 2 ) (A. 2) 

where cr is a representation dependent constant. A simpler parametrization of the phase 
space is to use 



dQ = sm6d6d(j) (A.3) 

with the normalizations J dQ = An, Y2 a Q a Q a = 1 and J dQQ-Q = An. The SU(2) Casimir 
is then restored by inspection. 

APPENDIX B: PROJECTION METHOD 

If we define the phase space density 5f™(kp,t) 

N 

5fr(kp, t)=Y, e- lkMt) S(p - Pl (t))Yi m m - nfo(p)S lo 5 m0 5 ko Y (B.l) 
we can construct structure factor Si(t, k,pp') for /th partial wave 



An -Ys( 6 fr(kp,t)\5fr(k P ',o)) 



21 + 



An 



i E 5>"* ( * ( * ) " r ' ( * )) *(p - Pi(t))s( P ' - p^rmrrm) - ^imhp') 
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m i, j 



= SKt,fe,Pp') (B.2) 
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Here a scalar product (A\B) is defined as (A*B) eq . We follow [15|, |22|, |23j and recast the 
formal Liouville equation flIII.4j) in the form of a formal eigenvalue-like equation in phase 
space 



S,(fcz;j>lO = (Sfr(kp)\(z - C)- X \8f?{k'p')) (B.3) 

The color charge effect by partial waves is represented as l,m in Eq. ( IB. 31) . If we introduce 
the projection operator 

V = 4nJ2 J dp 1 dp 2 \8fr(k,p 1 ))S^(k,p 1 ,p 2 )(8fr(k,p 2 )\ = 1 - Q (B.4) 

l,m,k 

we can check that this projection operator satisfies V 2 = V 



V 2 =AnJ2 E dp 1 dp 2 dp' 1 dp' 2 \8fr(k,p 1 ))S^(k,p 1 ,p 2 ) 

l,m,hl> , m ' ,k' 

x 47r(5/r(fe>P 2 )l^'(fe ) P / i))So/(fe,Pi',P2 / )(^'*(fe>P2 / )l = V (B.5) 
because of the translational invariance in space and the rotational invariance in color space, 



An(5fr(k )P2 )\5ff(k,p\)) = 5 kk/ 5 w 8 mm S ol (k, P2 , Pl ') (B.6) 

The off-diagonal elemenets vanish in the equilibrium averaging due to phase incoherence. 
Therefore, the projection operator in Eq. flB.5j) satisfies also Q 2 = Q and VQ = QV = 0. 
If we define |F z m (fep; z)) as |F z m (fep; 2)) = (z - C)- 1 ^ fp{kp)) from Eq. (|R3]), we have 



P(«-£)|if»(fcp;«)) =V\6fr(kp)) (B.7) 

P in Eq. (IB.5P is the operator which projects phase space function of a multipartle state 
with Z'th partial wave into a single particle state of the same parial wave, \5 fJJ 1 ' (kp)) , 
P\gf(kp)) = \5jf{kp)). Therefore Q^f^^kp)) = (1 - V)\6fr(kp)) = 0. With these 
in mind, we can modify the above equation further using V + Q = I 
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{Vz _ V cv - VCQ)\Fr(kp; z)) = V\6fr(kp)) 

(Qz - QCV - QCQ)\Fr(kp; z)) = (B.8) 
From these equations, we can extract 

zV\Fr(kp; z)) - VCV\Fr(kp; z)) - VCQ(z - QCQy 1 QCV\Fl n (kp; z)) = V\5fT(kp)) 

(B.9) 

By multiplying (5f(kp)\ we finally obtain, 

zSi(kz-pp') - J dp 1 dE l (kz;pp 1 )S l (kz;p 1 p') = S z (fcO;pp') (B.10) 
where the memory function, or the evolution operator Hi(kz; ppi) is 

Ukz; PP ') = / ^iWr(fc,P)l^ + ^l^ m (fc,Pi))S 7(fc,Pi,p') (B.H) 



with 



q? = VCQ(z - QCQ)~ 1 QCV (B.12) 



Since the Liouville operator C can be split into £ + £>i + £q, Eq. (111.9j) . the evolution 



operator can also be split into four terms; the free streaming term(£°), the self consistent 
term(Ef), the color charge term(S^) and the non-local collision term(E c ). 



k m p 

T, ol (kz;pp') = $(p-p') 

m 

k ' p 

E n (kz;pp f ) = -n fo(p)c D i(k) 

m 

Z Ql (kz;pp') = 

Y, C i{kz-pp>) = nf l M 2 ^ x j2(6fr(k P )\£Q(z - QCQr'Qmnkp')) (B.13) 
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APPENDIX C: COLLISIONAL COLOR CONTRIBUTION 



In this Appendix we detail the calculation that leads to a zero contribution from the 
colored Liouville operator in the collisional part of the self energy in the free streaming 
approximation. A typical contribution to (1V.2j) and (IV. 5 ft is 



Lg{q, «i) L%{q', q 2 ) $(t, q, q 2 ) S(t, q', q x ) = ^ 



x ( V(r -r 1 )QxQ 1 - (V Q - V Q 
x fc' D (r' - r 2 , Q' ■ Q 2 )) Q'xQ 2 - (V Q , - V Qa )) S(f, q, q 2 ) S(t, q', q x ) (C.l) 
which can be reduced to 

L Q (q, q x ) L*(q', q 2 ) S(t, q, q 2 ) S(t, q', q x ) = -- V(r - n) c' D {r' - r 2 , Q' ■ Q 2 j) 

x (s'(Q ■ Q 2 )S'(Q' • Q 1 )(Q 1 x Q 2 ) • Q (Q 1 x Q 2 ) ■ Q' 

xS"(Q ■ Q 2 )S(Q' ■ Q 1 )(Q 1 x Q 2 ) ■ Q (Q x Q') • Q 2 
xS(Q • Q 2 )S"{Q' ■ Q 1 )(Q x Q') • Q : (Q x x Q 2 ) ■ Q' 

xS'(Q • Q 2 )S'(Q' ■ Qx)(Q x Q') Q X {Q x Q') • Q 2 j (C.2) 

The derivatives onco and S are on their color argument. We note that (IC.2I) contribute to 
the collisional part of the self energy in (1B.6P after the integration over Q x and Q 2 , which 
is then zero. This is expected. Indeed, the colored Liouville operator is a 3-body force that 
requires 3 distinct color charges to not vanish. While (IC.2j) contributer to the unintegrated 
collisional operator, it does not in the integrated one which is the self-energy on the 2point 
function. It does contribute in the Liouville hierarchy in the 3-body structure factors and 
higher. 

APPENDIX D: HYDRO DYNAMICAL SUBSPACE 



The projection method onto the hydrodynamical subspace has been discussed by 



many 



id, 



11 



15] . This consists in dialing the projector in (I VI. 21) onto the hydrodynamical 
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modes. We choose Hermite polynomials as a basis set with the Maxwell-Boltzman distribu- 
tion fo(p) as a Gaussian weight function. The Hermite polynomials are the the generalized 
ones in 3D [25]. Specifically 

Hl(n) (p) = 1 #2(0 (p) = Pz #3(e) (p) = ^ (p 2 ~ 3) 

H 4(t 1 )(p)=Px H 5(t2 )(p)=Py (D.l) 

These polynomials are orthonormal for the inner product 
( m \n) = J dpa m H m (p)a n H n (p)nf (p) = S mn 

(m\F(k,t)\n) = J dpdp'a m H m (p)F(k,t;pp')a n H n (p')nf Q (p') (D.2) 

Here a m and a n set the normalizations. We chose the longitudinal momentum direction 
along k in Fourier space, (l\ = a m k ■ p. The transverse directional is chosen orthogonal to 
k, (t\ = a' m e ■ p with a unit vector satisfying e 2 = 1 and e ■ k = 0. 

The hydrodynamical projection operators Vh restricted to the five states fID.ip are 



5 5 

7^ = Qh = 1-V h = 1-Y,\W\ ( D - 3 ) 

i i 

While in general these 5 statesare enough to characterize the hydrodynamical modes in the 
SU(2) phase space, we need additional states to work out the shear viscosity as it involves 



261 ] . For that we need 



in general correlations in the stress tensor through the Kubo relation 
additionally, 



H 6 (p) = p x Py H 7 {p) = p x p z H 8 (p) = p y p z (D.4) 
With the definition of Gy(kz) = (i\S(kz; PP^inf^p))^^) we can rewrite flVI.2p as 



- Y J ^ kz ^PP')\ k )) G ^z) = Gy(feO) (D.5) 



k 



34 



where i,j are short for: n(density), e(energy), /(longitudinal momentum) and ^(transverse 
momentum). 
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